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Abstract. We explore the effect of zeros at the central point on nearby zeros of el- 
liptic curve L-functions, especially for one-pai'ameter families of rank r over Q. By the 
Birch and Swinnerton Dyer Conjecture and Silverman's Specialization Theorem, for t suf- 
ficiently large the L-function of each curve Et in the family has r zeros (called the family 
zeros) at the central point. We observe experimentally a repulsion of the zeros neai' the 
central point, and the repulsion increases with r. There is greater repulsion in the subset of 
curves of rank r + 2 than in the subset of curves of rank r in a rank r family. For curves 
with comparable conductors, the behavior of rank 2 curves in a rank one-parameter fam- 
ily over Q is statistically different from that of rank 2 curves from a rank 2 family. Unlike 
excess rank calculations, the repulsion decreases markedly as the conductors increase, and 
we conjecture that the r family zeros do not repel in the limit. Finally, the differences 
between adjacent normalized zeros near the central point are statistically independent of 
the repulsion, family rank and rank of the curves in the subset. Specifically, the differences 
between adjacent normalized zeros are statistically equal for all curves investigated with 
rank 0, 2 or 4 and comparable conductors from one-parameter families of rank or 2 over 

Q. 



1. Introduction 

Random matrix theory has successfully modeled the behavior of the zeros and values 
of many i-functions; see for example the excellent surveys IKeSn2l [Farl . The corre- 
spondence first appeared in Montgomery's analysis of the pair correlation' of the zeros of 
the Riemann zeta function as the zeros tend to infinity IMonl. Dyson noticed that Mont- 
gomery's answer, though limited to test functions satisfying certain support restrictions, 
agrees with the pair correlation of the eigenvalues from the Gaussian Unitary Ensemble^ 
(GUE). Montgomery conjectured that his result holds for all correlations and all support. 
Again with suitable restrictions and in the limit as the zeros tend to infinity, Hejhal |Hej| 
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'if {oLj is an increasing sequence of numbers and B C M""^ is a compact box, the n-level congelations 

are 

lim # 

JV-.00 TV 

One may replace the boxes with smooth test functions; see IRuSal for details. 

^The GUE is the N ^ oo scaling limit oi N X N complex Hermitian matrices with entries independently 
chosen from Gaussians; see IMehl for details. 
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showed the triple correlation of zeros of the Riemann zeta function agree with the GUE, 
and, more generally, Rudnick and Sarnak IRuSal showed the n-level correlations of the 
zeros of any principal L-function (the i-function attached to a cuspidal automorphic rep- 
resentation of GLm over Q) also agree with the GUE. 

In this paper we explore another connection between L-functions and random matrix 
theory, the effect of multiple zeros at the central point on nearby zeros of an L-function and 
the effect of multiple eigenvalues at 1 on nearby eigenvalues in a classical compact group. 
Particularly interesting cases are families of elliptic curve L-functions. It is conjectured 
that zeros of primitive L-functions are simple, except potentially at the central point for 
arithmetic reasons. For an elliptic curve E, the Birch and Swinnerton-Dyer Conjecture 
1BS-DlllBS^D2l states that the rank of the Mordell-Weil group E{Q) equals the order of 
vanishing of the L-function L{E, s) at the central point^ = 5- Let £ be a one-parameter 
family of elliptic curves over Q with (Mordell-Weil) rank'* r: 

= x^+A{T)x + B{T), A{T), B{T) e Z[T]. (1.1) 

For all t sufficiently large each curve Et in the family £ has rank at least r, by Silverman's 
specialization theorem ISi2l . Thus we expect each curve's L-function to have at least r 
zeros at the central point. We call the r conjectured zeros from the Birch and Swinnerton- 
Dyer Conjecture the family zeros. Thus, at least conjecturally, these families of elliptic 
curves offer an exciting and accessible laboratory where we can explore the effect of mul- 
tiple zeros on nearby zeros. 

The main tool for studying zeros near the central point (the low-lying zeros) in a family 
is the n-level density. Let (j){x) = Y[i=i 4>i{^i) where the (pi are even Schwartz functions 
whose Fourier transforms 0^ are compactly supported. Following Iwaniec-Luo-Sarnak 
IILSI . we define the n-level density for the zeros of an L-function L(.s, /) by 

D.A<t>) = E <t>Al,n'^)---^r.U,}^)-, (1.2) 

Cf is the analytic conductor of L{s,f), whose non-trivial zeros are i + 67/^. Under 
GRH, the non-trivial zeros all lie on the critical line !R(s) = i, and thus jjj G K. As 
(j>i is Schwartz, note that most of the contribution is from zeros near the central point. 
The analytic conductor of an L-function normalizes the non-trivial zeros of the L-function 
so that, near the central point, the average spacing between normalized zeros is 1; it is 
determined by analyzing the F-factors in the functional equation of the L-function (see for 
example IILSI V For elliptic curves the analytic conductor is the conductor of the elliptic 
curve (the level of the corresponding weight 2 cuspidal newform from the Modularity 
Theorem of IWillTaWillBCDTi V 

We order a family T of L-functions by analytic conductors. Let = {f ^ ^ '■ C f < 
iV}. The n-level density for the family T with test function is 

DnA^) - lim Dn,:FM). (1.3) 

where 

We normalize all L-functions to have functional equation s i— > 1 — s, and thus central point is at s = ^ ■ 
"^The group of rational function solutions (a:(T), y{T)) e iQ)(T)^ to y'^ = + A{T)x + B{T) is isomor- 
phic to Zi^ ® T, where T is the torsion part and r is the rank. 
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We can of course investigate other subsets. Other common choices are {/ : C/ £ 
[N,2N]}, or, for a one-parameter family £ of elliptic curves over Q, {Et E £ : t G 
[iV, 2A^]}. 

Let U(Af ) be the ensemble of N x N unitary matrices with Haar measure. The classical 
compact groups are sub-ensembles G{N) of U(A^); the most frequently encountered ones 
are USp(2Af), S0(2iV) and S0(2iV + 1). Katz and Sarnak's Density Conjecture IKaSall 
IKaSa2l states that as the conductors tend to infinity, the behavior of the normalized zeros 
near the central point equals the N ^ oo scaling limit of the normalized eigenvalues near 
1 of a classical compact group; see il.H for an exact statement. In the function field case, 
the corresponding classical compact group can be identified from the monodromy group; 
in the number field case, however, the reason behind the identification is often a mystery 
(see IDMI l As the eigenvalues of a unitary matrix are of the form e**, we often talk about 
the eigenangles 9 instead of the eigenvalues e'^, and the eigenangle corresponds to the 
eigenvalue 1. 

Using the explicit formula we replace the sums over zeros in (II. 2> with sums over the 
Fourier coefficients at prime powers. For example, if E : = + Ax + B is an elliptic 
curve, assuming GRH the non-trivial zeros of the associated i-function 

oc 

L{E,s) = Y.^E{n)n-' (1.5) 

n=l 

(normalized to have functional equation s i-^ 1 — s) are i + 17, 7 e K. If (/) is a Schwartz 
test function, then the explicit formula for L{E, s) is 



^U'^) = ^(0) + 0(0) -2^ .1°^^ .HogMA.b) 



V V 2^ / v^w rv ; ^ \ogCE \logCE J VP 



2E 



logp 21ogp\ Xj^{p) 



log Ce \log Ce 

see for example IMesllMilTl . By using appropriate averaging formulas and combinatorics, 
the resulting prime power sums in the n-level densities can be evaluated for of suitably 
restricted support. The Density Conjecture is that to each family of i-functions for any 
Schwartz test function 4> ■ K" ^ K", 

Dn.A(^) - j 4>{x)W„^g{x)dx = J ^{u)Wn,g{u)du. (1.7) 

The density kernel Wn.g{x) is determined from the N oo scaling limit of the asso- 
ciated classical compact group C{N); the last equality follows by Plancherel. The most 
frequently occurring answers are the scahng limits of Unitary, Symplectic and Orthogonal 
ensembles. For n = 1 we have 

Wi^viu) = S{u) 

m,vsp{u) = 5{u) - il(u) 

2l,SO(eve„)(«) = S{u) + ^l{u) (1.8) 

W^i,so(odd)(w) 5{u)~^I{u) + l 
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where l{u) is the characteristic function of [—1, 1]. For arbitrarily small support, unitary 
and symplectic are distinguishable from each other and the orthogonal groups; however, 
for test functions cfi supported in (—1, 1), the three orthogonal groups agree: 

/ (j){u)Wi^\j{u)du = <f){u) 

J ^{u)W,,usp{u)du = 0(u)-i(A(O) 

J liu)W,^SOie.en)iu)du = + i</,(0) (1.9) 

J Hu)Wi^so{odd)iu)du = 0(ii) + i0(O) 
J $iu)Wi,oiu)du = + i(/.(0). 

Similar results hold for the ri-level densities, but below we only need the 1-level; see 
IConllKaSall for the derivations of the general n-level densities, and Appendix|X|for the 
1-level density for the orthogonal groups. 

For one-parameter families of elliptic curves, the results suggest that the correct models 
are orthogonal groups (if all functional equations are even then the answer is SO(cvcn), 
while if all are odd the answer is SO (odd)). Often instead of normalizing each curve's ze- 
ros by the logarithm of its conductor (the local rescaling), one instead uses the average log- 
conductor (the global rescaling). If we are only interested in the average rank, it suffices to 
study just the 1-level density from the global rescaling. This is because we only care about 
the imaginary parts of the zeros at the central point, and both scalings of the imaginary 
part of the central point are zero; see for example IBrI IGo2l IH-BI iMicl ISi3l IYo2l . To date 
all results have support in (—1, 1), where il.9i shows that the behavior of O, SO(even) 
and SO (odd) are indistinguishable. If we want to specify a unique corresponding classical 
compact group we study the 2-level density as well, which for arbitrarily small support 
suffices to distinguish the three orthogonal candidates. Using the global rescaling removes 
many complications in the 1-level sums but not in the 2-level sums. In fact, for the 2- 
level investigations the global rescaling is as difficult as the local rescaling; see IMil2l for 
details. 

Our research was motivated by investigations on the distribution of rank in families of 
elliptic curves as the conductors grow. As we see below, for the ranges of conductors stud- 
ied there is poor agreement between elliptic curve rank data and the N ^ oo scaling limits 
of random matrix theory. The purpose of this research is to show that another statistic, the 
distribution of the first few zeros above the central point, converges more rapidly. 

We briefly review the excess rank phenomenon. A generic one-parameter family of 
elliptic curves over Q has half of its functional equations even and half odd (see IHel for 
the precise conditions for a family to be generic). Consider such a one-parameter family of 
elliptic curves over Q, of rank r, and assume the Birch and Swinnerton-Dyer Conjecture. 
It is believed that the behavior of the non-family zeros is modeled by the N ^ oo scaling 
limit of orthogonal matrices. Thus if the Density Conjecture is correct, then at the central 
point in the limit as the conductors tend to infinity the L-functions have exactly r zeros 
50% of the time, and exactly r + 1 zeros 50% of the time. Thus in the limit half the curves 
have rank r and half have rank r + 1. In a variety of families, however, one observes^ that 



Actually, this is not quite true. The analytic rank is estimated by the location of the first non-zero term in the 
series expansion of L(E, s) at the central point (see ICrl for the algorithms). For example, if the zeroth through 
third coefficients are smaller than 10~® and the fourth is 1.701, then we say the curve has analytic rank 4, even 
though it is possible (though unlikely) that one of the first four coefficients is really non-zero. It is difficult to 
prove an elliptic curve L-function vanishes to order two or greater. Goldfeld IGoll and Gross-Zagier IGZI give 
an effective lower bound for the class number of imaginary quadratic fields by an analysis of an elliptic curve 
L-function which is proven to have three zeros at the central point. 
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30% to 40% have rank r, about 48% have rank r + 1, 10% to 20% have rank r + 2, and 
about 2% have rank r + 3; see for example l5MllFeTllFe2IIZKI . 

We give a representative family below; see in particular IFe2l for more examples. Con- 
sider the one-parameter family y'^ — + 16Tx + 32 of rank over Q. Each range below 
has 2000 curves: 

T-range rank rank 1 rank 2 rank 3 run time (hours) 



[-1000, 1000) 39.4% 47.8% 12.3% 0.6% <1 

[1000, 3000) 38.4% 47.3% 13.6% 0.6% <1 

[4000, 6000) 37.4% 47.8% 13.7% 1.1% 1 

[8000, 10000) 37.3% 48.8% 12.9% 1.0% 2.5 

[24000,26000) 35.1% 50.1% 13.9% 0.8% 6.8 

[50000,52000) 36.7% 48.3% 13.8% 1.2% 51.8 



The relative stability of the percentage of curves in a family with rank 2 above the family 
rank r naturally leads to the question as to whether or not this persists in the limit; it cannot 
persist if the Density Conjecture (with orthogonal groups) is true for all support^. Recently 
Watkins IWatI investigated the family x^+y^ = m for varying m, and unlike other families 
his range of m was large enough to see the percentage with rank r + 2 markedly decrease, 
providing support for the Density Conjecture (with orthogonal groups). 

In our example above, as well as the other families investigated, the logarithms of the 
conductors are quite small. Even in our last set the log-conductors are under 40. An analy- 
sis of the error terms in the explicit formula suggests the rate of convergence of quantities 
related to zeros of elliptic curves is like the logarithm of the conductors. It is quite sat- 
isfying when we study the first few normalized zeros above the central point that, unlike 
excess rank, we see a dramatic decrease in repulsion with modest increases in conductor 

In ^we study two random matrix ensembles which are natural candidates to model 
families of elliptic curves with positive rank. Many natural questions concerning the nor- 
malized eigenvalues for these models for finite N lead to quantities that are expressed in 
terms of eigenvalues of integral equations. Our hope is that showing the possible connec- 
tions between these models and number theory will spur interest in studying these models 
and analyzing these integral equations. We assume the Birch and Swinnerton-Dyer Con- 
jecture, as well as GRH. We calculate some properties of these ensembles in Appendix 

m 

In ^we summarize the theoretical results of previous investigations, which state: 

• For one-parameter families of rank r over Q and suitably restricted test functions, 
as the conductors tend to infinity the 1 -level densities imply that in this restricted 
range, the r family zeros at the central point are independent of the remaining 
zeros. 

If this were to hold for all test functions, then as the conductors tend to infinity the 
distribution of the first few normalized zeros above the central point would be independent 
of the r family zeros. 

In 13 we numerically investigate the first few normalized zeros above the central point 
for elliptic curves from many families of different rank. Our main observations are: 

• The first few normalized zeros are repelled from the central point. The repulsion 
increases with the number of zeros at the central point, and even in the case when 



Explicitly, if the large-conductor limit of the elliptic curve L-functions agree with the — > oo scaling limits 
of orthogonal groups. 
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there are no zeros at the central point there is repulsion from the large-conductor 
limit theoretical prediction. This is observed for the family of all elliptic curves, 
and for one-parameter families of rank r over Q. 

• There is greater repulsion in the first normalized zero above the central point for 
subsets of curves of rank 2 from one-parameter families of rank over Q than 
for subsets of curves of rank 2 from one-parameter families of rank 2 over Q. It 
is conjectured that as the conductors tend to infinity, 0% of curves in a family of 
rank r have rank r + 2 or greater. If this is true, we are comparing a subset of zero 
relative measure to one of positive measure. As the first set is (conjecturally) so 
small, it is not surprising that to date there is no known theoretical agreement with 
any random matrix model for this case. 

• Unlike most excess rank investigations, as the conductors increase the repulsion 
of the first few normalized zeros markedly decreases. This supports the conjecture 
that, in the limit as the conductors tend to infinity, the family zeros are indepen- 
dent of the remaining normalized zeros (i.e., the repulsion from the family zeros 
vanishes in the limit). 

• The repulsion from additional zeros at the central point cannot entirely be ex- 
plained by collapsing some zeros to the central point and leaving all the other 
zeros alone. See in particular Remark l431 

• While the first few normalized zeros are repelled from the central point, the differ- 
ences between normalized zeros near the central point are statistically independent 
of the repulsion, as well as the method of construction. Specifically, the differences 
between adjacent zeros near the central point from curves of rank 0, 2 or 4 with 
comparable conductors from one-parameter families of rank or 2 over Q are sta- 
tistically equal. Thus the data suggests that the effect of the repulsion is simply to 
shift aU zeros by approximately the same amount. 

The numerical data is similar to excess rank investigations. While both seem to con- 
tradict the Density Conjecture, the Density Conjecture describes the limiting behavior as 
the conductors tend to infinity. The rate of convergence is expected to be on the order of 
the logarithms of the conductors, which is under 40 for our curves. Thus our experimental 
results are likely misleading as to the limiting behavior It is quite interesting that, unlike 
most excess rank investigations, we can easily go far enough to see conductor dependent 
behavior. 

Thus our theoretical and numerical results, as well as the Birch and Swinnerton-Dyer 
and Density Conjectures, lead us to 

Conjecture 1.1. Consider one-parameter families of elliptic curves of rank r over Q and 
their sub-families of curves with rank exactly r + k for k G {0, 1,2,...}. For each sub- 
family there are r family zeros at the central point, and these zeros repel the nearby nor- 
malized zeros. The repulsion increases with r and decreases to zero as the conductors tend 
to infinity, implying that in the limit the r family zeros are independent of the remaining 
zeros. Ifk>2 these additional non-family zeros at the central point may influence nearby 
zeros, even in the limit as the conductors tend to infinity. The spacings between adjacent 
normalized zeros above the central point are independent of the repulsion; in particular, it 
does not depend on r or k, but only on the conductors. 

2. Random Matrix Models for Families of Elliptic Curves 

We want a random matrix model for the behavior of zeros from families of elliptic 
curve i-functions with a prescribed number of zeros at the central point. We concentrate 
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on models for either one or two-parameter families over Q, and refer the reader to IIFarl for 
more on random matrix models. Both of these families are expected to have orthogonal 
symmetries. Many people (see for example IIDFKI IUo2l IGMI IMa i' 'RuSi' ' Rub2l ISTI ) have 
studied famiUes constructed by twisting a fixed elliptic curve by characters. The general 
belief is that such twisting should lead to unitary or symplectic families, depending on the 
orders of the characters. 

There are two natural models for the corresponding situation in random matrix theory 
of a prescribed number of eigenvalues at 1 in sub-ensembles of orthogonal groups. For 
ease of presentation we consider the case of an even number of eigenvalues at 1; the odd 
case is handled similarly. 

Consider a matrix in S0{2N). It has 2N eigenvalues in pairs e^^^' , with 6j G [0, n]. 
The joint probability measure on e = (01,...,0Ar) e [0,7r]^is 

deo(e) = CAT Y[ (cose'fc - cos6lj)2 J| dOj, (2.10) 

l<j<k<N i<j<N 

where cjv is a normalization constant so that deo{'d) integrates to 1. From ( I2.10t we can 
derive all quantities of interest on the random matrix side; in particular, n-level densities, 
distribution of first normalized eigenvalue above 1 (or eigenangle above 0), and so forth. 

We now consider two models for sub-ensembles of S0{2N) with 2r eigenvalues at 1, 
and the N oo scaling limit of each. 

Independent Model: The sub-ensemble of S0(2iV) with the upper left block a 2r x 2r 
identity matrix. The joint probability density of the remaining N — r pairs is given by 

de2r,Indop(B) = C2rJndcp.Ar ]J (cOS 6*^ - COS )^ ]J dO j . (2.11) 

l<j<fe<7V-r l<i<Ar-r 

Thus the ensemble is matrices of the form 

^ :.geSO(2iV-2r)|; (2.12) 

the probabilities are equivalent to choosing g with respect to Haar measure on S0{2N — 
2r). We call this the Independent Model as the forced eigenvalues at 1 from the I2rx2r 
block do not interact with the eigenvalues of g. In particular, the distribution of the re- 
maining N — r pairs of eigenvalues is exactly that of S0(2iV — 2r); this block's N ^ oo 
scaling limit is just SO(cvcn). See IConllKaSall as well as Appendix IaI 

Interaction Model: The sub-ensemble of S0(2iV) with 2r of the 27V eigenvalues 
equaling 1: 

de2rJntor(B) = C2rJntGr,JV J| {cOS 6 k - COS 6 j)'^ J]^ {1 - COS 9 j)'^'' dO j . 

l<j<k<N-r l<j<N-r 

(2.13) 

We call this the Interaction Model as the forced eigenvalues at 1 do affect the behavior of 
the other eigenvalues near 1. Note here we condition on all S0{2N) matrices with at least 
2r eigenvalues equal to 1. The (1 — cos 9j)'^'^ factor results in the forced eigenvalues at 1 
repelling the nearby eigenvalues. 

Remark 2.1. As the calculations for the local statistics near the eigenvalue at 1 in the 
Interaction Model has not appeared in print, in AppendixlAl (written by Eduardo Duenez) 
is a derivation of formula ( I2.13> (see especially ^A.2t . as well as the relevant integral 
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(Bessel) kernels dictating such statistics. See also ISnl for the value distribution of the first 
non-zero derivative of the characteristic polynomials of this ensemble. 

While both models have at least 2r eigenvalues equal to 1, they are very different sub- 
ensembles of S0{2N), and they have distinct limiting behavior (see also Remark im . 
We can see this by computing the 1-level density for each, and comparing with ( ll.9l l. 

Letting Tyi,so(ovcn) (respectively Vl^i,so(cvcn),indcp,2r and W'i,so(ovon),intcr,2r) denote the 
Fourier transform of the kernel for the 1-level density of SO(even) (respectively, of the 
Independent Model for the sub-ensemble of SO (even) with 2r eigenvalues at 1 and of the 
Interaction Model for the sub-ensemble of SO (even) with 2r eigenvalues at 1), we find in 
AppendixlAlthat 

W^l,SO(ovon)(w) = S{u) + ^l{u) 
M^l,SO(ovcn),Indcp,2r(w) = S{u) + ^l{u) + 2 

W'l,SO(cvc„)Jntcr,2r(u) = S (u) + ^l{u) + 2 + 2 {\u\ - l)l{u) . (2.14) 

As l{u) is positive for \u\ < 1, note that the density is smaller for \u\ < 1 in the Interaction 
versus the Independent Model. We can interpret this as a repulsion of zeros, as the follow- 
ing heuristic shows (though see Appendix IaI for proofs). We compare the 1-level density 
of zeros from curves with and without repulsion, and show that for a positive decreasing 
test function, the 1-level density is smaller when there is repulsion. 

Consider two elliptic curves, E of rank and conductor Ce and E' of rank r and 
conductor Ce'- Assume Ce ^ Ce' ^ C, and assume GRH for both i-functions. If the 
curve E has rank then we expect the first zero above the central point, ^ + 67^; i, to have 
^E,i ~ logC ' ^'^^ ^ family zeros at the central point repel, it is reasonable to 

posit a repulsion of size -^^q for some br > 0. This is because the natural scale for the 
distance between the low-lying zeros is i^^q , so we are merely positing that the repulsion 
is proportional to the distance. We assume all zeros are repelled equally; evidence for this 
is provided in ^4.61 Thus for E' (the repulsion case) we assume "/e'.j ~ lE.j + k^^- 
We can detect this repulsion by comparing the 1-level densities of E and E' . Take a non- 
negative decreasing Schwartz test function 4>. The difference between the contribution 
from the j* zero of each is 

logC\ ,/ logC\ ,/ logC br\ ,/ logC 



^ ( \ogC\ br 

" ^(^^^..— j-^- (2.15) 

As (p is decreasing, its derivative is negative and thus the above shows the 1-level density 
for the zeros from E' (assuming repulsion) is smaller than the 1 -level density for zeros from 
E. Thus the lower 1-level density in the Interaction Model versus the Independent Model 
can be interpreted as a repulsion; however, this repulsion can be shared among several 
zeros near the central point. In fact, the observations in ^4.61 suggest that the repulsion 
shifts all normalized zeros near the central point approximately equally. 

3. Theoretical Results 



Consider a one-parameter family of elliptic curves of rank r over Q. We summarize 
previous investigations of the effect of the (conjectured) r family zeros on the other zeros 
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near the central point. For convenience we state the results for the global rescaling, though 
similar results hold for the local rescaUng (under slightly more restrictive conditions; see 
IIMil2l for details). For small support, the 1 and 2-level densities agree with the scaling 
limits of 

(Iry.r \ j Irxr \ I Irxr \ ("I 1 

0{N) J' \ S0{2N) J' y S0(2iV + l) J' ^^'^^^ 

depending on whether or not the signs of the functional equation are equidistributed or all 
the signs are even or all the signs are odd. The 1 and 2-level densities provide evidence 
towards the Katz-Sarnak Density Conjecture for test functions whose Fourier transforms 
have small support (the support is computable and depends on the family). See IMilll 
for the calculations with the global rescaling, though the result for the 1-level density is 
implicit in IISi3l . Similar results are observed for two-parameter families of elliptic curves 
in |'MiTr Yo2l. 

While the above results are consistent with the Birch and Swinnerton-Dyer Conjecture 
that each curve's L-function has at least r zeros at the central point, it is not a proof (even 
in the limit) because our supports are finite. For families with t £ [N, 2N] the errors are 
of size Ql iogjy ) or 0{ ^°f^°fj^ ). Thus for large N we cannot distinguish a family with 
exactly r zeros at the central point from a family where each Et has exactly r zeros at 

For one-parameter families of elliptic curves over Q, in the limit as the conductors tend 
to infinity the family zeros (those arising from our belief in the Birch and Swinnerton-Dyer 
Conjecture) appear to be independent from the other zeros. Equivalently, if we remove the 
contributions from the r family zeros, for test functions with suitably restricted support 
the spacing statistics of the remaining zeros agree perfectly with the standard orthogonal 
groups O, SO (even) and SO (odd), and it is conjectured that these results should hold for 
all support. Thus the n-level density arguments support the Independent over the Inter- 
action Model when we study all curves in a family; however, these theoretical arguments 
do not apply if we study the sub-family of curves of rank r + k {k > 2) in a rank r 
one-parameter family over Q. 

Remark 3.1. It is important to note that our theoretical results are for the entire one- 
parameter family. Specifically, consider the subset of curves of rank r + 2 from a one- 
parameter family of rank r over Q. If the Density Conjecture (with orthogonal groups) is 
true, then in the limit 0% of curves are in this sub-family. Thus these curves may behave 
differently without contradicting the theoretical results for the entire family. Situations 
where sub-ensembles behave differently than the entire ensemble are well known in ran- 
dom matrix theory. For example, to any simple graph we may attach a real symmetric 
matrix, its adjacency matrix, where aij = 1 if there is an edge connecting vertices i and 
j, and otherwise. The adjacency matrices of d-regular graphs are a thin sub-ensemble of 
real symmetric matrices with entries independently chosen from { — 1, 0, 1}. The density 
of normalized eigenvalues in the two cases are quite different, given by Kesten's Measure 
IMcKI for d-regular graphs and Wigner's Semi-Circle Law IMehl for the real symmetric 
matrices. 

It is an interesting question to determine the appropriate random matrix model for rank 
r + 2 curves in a rank r one-parameter family over Q, both in the limit of large conductors 
as well as for finite conductors. We explore this issue in greater detail in ^4.31 to ^4.61 
where we compare the behavior of rank 2 curves from rank one-parameter families over 
Q to that of rank 2 curves from rank 2 one-parameter families over Q. 
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4. Experimental Results 

We investigate the first few normalized zeros above the central point. We used Michael 
Rubinstein's L-function calculator IRub3l to determine the zeros. The program does a 
contour integral to ensure that we found all the zeros in a region, which is essential in 
studies of the first zero! See IRubll for a description of the algorithms. The analytic 
ranks were found (see Footnote [Q by determining the values of the L-functions and their 
derivatives at the central point by the standard series expansion; see ICrl for the algorithms. 
Some of the programs and all of the data (minimal model, conductor, discriminant, sign of 
the functional equation, first non-zero Taylor coefficient from the series expansion at the 
central point, and the first three zeros above the central point) are available onhne at 

http : / / www . math . brown .edu/^sj miller/ repulsion 

We study several one-parameter families of elliptic curves over Q. As all of our fami- 
lies are rational surfaces^, Rosen and Silverman's result that the weighted average of fibral 
Frobenius trace values determines the rank over Q (see IRoSil ) is applicable, and evalu- 
ating simple Legendre sums suffices to determine the rank. We mostly use one-parameter 
families from Fermigier's tables IIFe2l . though see lALMI for how to use the results of 
IIRoSil to construct additional one-parameter families with rank over Q. 

We cannot obtain a decent number of curves with approximately equal log-conductors 
by considering a solitary one-parameter family. The conductors in a family typically grow 
polynomially in t. The number of Fourier coefficients needed to study a value of L{s,Et) 
on the critical line is of order \/Ct\ogCt (Ct is the conductor of Et), and we must then 
additionally evaluate numerous special functions. We can readily calculate the needed 
quantities up to conductors of size 10^^, which usually translates to just a few curves in a 
family. We first studied all elliptic curves (parametrized with more than one parameter), 
found the minimal models, and then sorted by conductor We then studied several one- 
parameter families, amalgamating data from different families if the curves had the same 
rank and similar log-conductor. 

Remark 4.1. Amalgamating data from different one-parameter families warrants some 
discussion. We expect that the behavior of zeros from curves with similar conductors and 
the same number of zeros and family zeros at the central point should be approximately 
equal. In other words, we hope that curves with the same rank and approximately equal 
conductors from different one-parameter families of the same rank r over Q behave simi- 
larly, and we may regard the different one-parameter families of rank r over Q as different 
measurements of this universal behavior. This is similar to numerical investigations of the 
spacings of energy levels of heavy nuclei; see for example IHHI IHP51 . In studying the 
spacings of these energy levels, there were very few (typically between 10 and 100) levels 
for each nucleus. The belief is that nuclei with the same angular momentum and parity 
should behave similarly. The resulting amalgamations often have thousands of spacings 
and excellent agreement with random matrix predictions. 

Similar to the excess rank phenomenon, we found disagreement between the experi- 
ments and the predicted large-conductor limit; however, we believe that this disagreement 
is due to the fact that the logarithms of the conductors investigated are small. In ^4.21 to 
^4.5l we find that for curves with zeros at the central point, the first normalized zero above 



''An elliptic surface y'^ = + A{T)x + B{T), A{T), B{T) G Z[T], is a rational surface if and only 
if one of the following is true: (1) < max{3degA, 2degB} < 12; (2) 3degA = 2deg_B = 12 and 
ordT=ori^A(T-i) = 0. 
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the central point is repelled, and the more zeros at the central point, the greater the re- 
pulsion. However, the repulsion decreases as the conductors increase. Thus the repulsion 
is probably due to the small conductors, and in the limit the Independent Model (which 
agrees with the function field analogue and the theoretical results of ^ should correctly 
describe the first normalized zero above the central point in curves of rank r in families of 
rank r over Q. It is not known what the correct model is for curves of rank r + 2 in a fam- 
ily of rank r over Q, though it is reasonable to conjecture it is the Interaction model with 
the sizes of the matrices related to the logarithms of the conductors. Keating and Snaith 
IKeSnlllKeSn2l showed that to study zeros at height T it is better to look at N x N matri- 
ces, with N = logT, than to look at the — > oo scaling limit. A fascinating question is 
to determine the correct finite conductor analogue for the two different cases here. Inter- 
estingly, we see in ^4.6l that the spacings between adjacent normalized zeros is statistically 
independent of the repulsion, which implies that the effect of the zeros at the central point 
(for finite conductors) is to shift all the nearby zeros approximately equally. 

4. 1 . Theoretical Predictions: Independent Model. In Figures 1 and 2 we plot the first 
normalized eigenangle above for SO(27V) (i.e., SO(even)) and S0(2iV+l) (i.e., SO(odd)) 
matrices. The eigenvalues occur in pairs e^^^' , 9j e [0, tt] ; by normalized eigenangles for 
SO (even) or SO (odd) we mean 6^^. We chose 2A^ < 6 and 2A^ + 1 = 7 for our 
simulations, and chose our matrices with respect to the appropriate Haar measure*^. We 
thank Michael Rubinstein for sharing his — > cx) scaling limit plots for S0(2A^) and 
S0(2iV + 1). 



Note that for SO(odd) matrices there is always an eigenvalue at 1. The N —> oo scaling limit of the 
distribution of the second eigenangle for SO (odd) matrices equals the A' — > oo scaling limit of the distribution 
of the first eigenangle for USp (Unitary Symplectic) matiices; see pages 10-11 and 411^16 of IKaSall and 
page 10 of IKaSa2l . 
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Figure la: First normalized eigenangle above 0: 23,040 SO(4) matrices 
Mean = .357, Standard Deviation about the Mean = .302, Median = .357 




Figure lb: First normalized eigenangle above 0: 23,040 SO(6) matrices 
Mean = .325, Standard Deviation about the Mean = .284, Median = .325 




Figure Ic: First normalized eigenangle above 0: 
N ^ oo scaling Umit of SO(2A^): Mean = .321. 
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Figure 2a; First normalized eigenangle above 1: 322,560 SO(7) matrices 
Mean = .879, Standard Deviation about the Mean = .361, Median = .879 




o75 i TTs 2 2T5 

Figure 2b: First normalized eigenangle above 1: 
^ 00 scaling limit of SO(2A^ + 1): Mean = .782. 

For the S0(2A^) matrices, note the mean decreases as 2N increases. A similar result 
holds for S0(2A^ + 1) matrices; as we primarily study even rank below, we concentrate on 
S0(2iV) here. As TV ^ 00, Katz and Sarnak (pages 412-415 of IIKaSa2l ) prove that the 
mean of the first normalized eigenangle above 6 — (corresponding to the eigenvalue 1) 
for SO(cven) is approximately 0.321, while for SO(odd) it is approximately 0.782. 

We study the first normalized zero above the central point for elliptic curve L-functions 
in ^.2l to ^4.51 We rescale each zero: '-)Et.i ^ lEt.i '"fif* • mean of the first normal- 
ized eigenangle above for S0(2A^) matrices decreases as 2N increases, and similarly 
we see that the first normalized zero above the central point in families of elHptic curves 
decreases as the conductor increases. This suggests that a good finite conductor model for 
families of elliptic curves with even functional equation and conductors of size C would 
be S0(2A^), with N some function of log C. 
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4.2. All Curves. 

4.2.1. Rank Curves. We study the first normalized zero above the central point for 
1500 rank elliptic curves, 750 with log(cond) e [3.2, 12.6] in Figure 3 and 750 with 
log(cond) e [12.6, 14.9] in Figure 4. These curves were obtained as follows: an elliptic 
curve can be written in Weierstrass form as 

y'^ + aixy + a^y = + a2x'^ + a4X + , a.; G Z. (4.17) 

We often denote the curve by [ai, a2, 03, 04, ag]. We let ai range from to 10 (as without 
loss of generality we may assume ai > 0) and the other a; range from —10 to 10. We 
kept only non-singular curves. We took minimal Weierstrass models for the ones left, and 
pruned the list to ensure that all the remaining curves were distinct. We then analyzed the 
first few zeros above the central point for 1500 of these curves (due to the length of time it 
takes to compute zeros for the curves, it was impossible to analyze the entire set). 

Figures 3 and 4 suggest that as the conductor increases the repulsion decreases. For the 
larger conductors in Figure 4, the results are closer to the predictions of Katz-Sarnak, and 
the shape of the distribution with larger conductors is closer to the random matrix theory 
plots of Figure 1. Though both plots in Figure 3 and 4 differ from the random matrix the- 
ory plots, the plot in Figure 4 is more peaked, the peak occurs earlier, and the decay in the 
tail is faster Standard statistical tests show the two means (1.04 for the smaller conduc- 
tors and 0.88 for the larger) are significantly different. Two possible tests are the Pooled 
Two-Sample t-Procedure"^ (where we assume the data are independently drawn from two 
normal distributions with the same mean and variance) and the Unpooled Two-Sample 
<-Procedure'" (where we assume the data are independently drawn from two normal distri- 
butions with the same mean and no assumption is made on the variance). See for example 
ICaBel . pages 409-410. Both tests give t-statistics around 10.5 with over 1400 degrees of 
freedom. As the number of degrees of freedom is so large, we may use the Central Limit 
Theorem and replace the t-statistic with a z-statistic. As for the standard normal the proba- 
bility of being at least 10.5 standard deviations from zero is less than 3.2 x 10^^^ percent, 
we obtain strong evidence against the null hypothesis that the two means are equal (i.e., 
we obtain evidence that the repulsion decreases as the conductor increases). 



^The Pooled Two-Sample t-Procedure is 



t = {Xi-X2)/sp.l — + —, (4.18) 
where Xi is the sample mean of observations of population i, Si is the sample standard deviation and 



(ni - l)s? + (n2 - l)s? 

^^J^ ' 1^^ (4 19) 

"1 + "2 - 2 



is the pooled variance; t has a t-distribution with rii + n2 — 2 degrees of freedom. 
''^Notation as in Footnotel9l the Unpooled Two-Sample t-Procedure is 



^2 j,2 

t = (Xi-X2)/i/-^ + ^; (4.20) 
ni 722 



this is approximately a t distribution with 

(m - 1) (n2 - 1) {n2sl +nis|)^ 
(n2 — l)'^2''l + ~ ^)"'1^2 

degrees of freedom 



(4.21) 
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Figure 3: First normalized zero above the central point: 
750 rank curves from + ai_xy + a-^y = + a2X^ + a^x 
log(cond) e [3.2, 12.6], median = 1.00, mean = 1.04, 
standard deviation about the mean = .32 
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Figure 4: First normalized zero above the central point: 
750 rank curves from j/^ + aixy + aj,y = x^ + a2X^ + aiX + ae, 
log(cond) e [12.6, 14.9], median = .85, mean = .88, 
standard deviation about the mean = .27 



4.2.2. Rank 2 Curves. We study the first normalized zero above the central point for 1330 
rank 2 elliptic curves, 665 with log(cond) e [10, 10.3125] in Figure 5 and 665 with 
log(cond) e [16, 16.5] in Figure 6. These curves were obtained from the same proce- 
dure which generated the 1500 curves in ^.2.11 except now we chose 1330 curves with 
what we believe is analytic rank exactly 2. We did this by showing the L-function has even 
sign, the value at the central point is zero to at least 5 digits, and the second derivative at 
the central point is non-zero; see also Footnote[2 In ^4.3l and ^4.4l we study other families 
of curves of rank 2 (rank 2 curves from rank and rank 2 one-parameter families over Q). 

The results are very noticeable. The first normalized zero is significantly higher here 
than for the rank curves. This supports the belief that, for small conductors, the repulsion 
of the first normalized zero increases with the number of zeros at the central point. 

We again split the data into two sets (Figures 5 and 6) based on the size of the conductor 
As the conductors increase the mean (and hence the repulsion) significantly decreases, 
from 2.30 to 1.82. 

We are investigating rank 2 curves from the family of all elliptic curves (which is a 
many parameter rank family). In the limit we believe half of the curves are rank and 
half are rank 1. The natural question is to determine the appropriate model for this subset 
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of curves. As in the limit we believe a curve has rank 2 (or more) with probability 
this is a question about conditional probabilities. 




1.5 2 2.5 3 3.5 

Figure 5: First normalized zero above the central point: 
665 rank 2 curves from + a\xy + a-^y = + 02^^ + a^x + ag. 
log(coiid) e [10, 10.3125], median = 2.29, mean = 2.30 
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Figure 6: First normalized zero above the central point: 
665 rank 2 curves from + aixy + a^y = + a2X^ + a^x 
log(cond) e [16, 16.5], median — 1.81, mean = 1.82 
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4.3. One-Parameter Families of Rank Over Q. 

4.3. 1 . Rank Curves. We analyzed 14 one-parameter families of rank over Q; we chose 
these families from IIFe2l . We want to study rank curves in a solitary one-parameter 
family; however, the conductors grow rapidly and we can only analyze the first few zeros 
from a small number of curves in a family. For our conductor ranges it takes several hours 
of computer time to find the first few zeros for all the curves in a family. In Figures 7 and 
8 and Tables 1 and 2 we study the first normalized zero above the central point for 14 one- 
parameter families of elliptic curves of rank over Q. Even though we have few data points 
in each family, we note the medians and means are always higher for the smaller conductors 
than the larger ones. Thus the "repulsion" is decreasing with increasing conductor, though 
perhaps repulsion is the wrong word here as there is no zero at the central point! We 
studied the median as well as the mean because, for small data sets, one or two outliers can 
significantly affect the mean; the median is more robust. 

For both the Pooled and Unpooled Two-Sample t-Procedure the ^-statistic exceeds 20 
with over 200 degrees of freedom. The Central Limit Theorem is an excellent approxi- 
mation and yields a z-statistic exceeding 20, which strongly argues for rejecting the null 
hypothesis that the two means are equal (i.e., providing evidence that the repulsion de- 
creases with increasing conductors). Note the first normalized zero above the central point 
is significantly larger than the N oo scaling limit of S0(2A^) matrices, which is about 
0.321. 

Some justification is required for regarding the data from the 14 families as independent 
samples from the same distribution. It is possible that there are family-specific lower order 
terms to the n-level densities (see LMil 1 Mil3 Yo2 j). Our amalgamation of the data is 
similar to physicists combining the energy level data from different heavy nuclei with 
similar quantum numbers. The hope is that the systems are similar enough to justify such 
averaging as it is impractical to obtain sufficient data for just one nucleus (or one family of 
elliptic curves, as we see in ^4.4> . See also Remark lTTI 

Remark 4.2. The families are not independent: there are 11 curves that occur twice and 
one that occurs three times in the small conductor set of 220 curves, and 133 repeats in 
the large conductor set of 996 curves. In our amalgamations of the families, we present 
the results when we double count these curves as well as when we keep only one curve in 
each repeated set. In both cases the repeats account for a sizeable percentage of the total 
number of observations; however, there is no significant difference between the two sets. 
Any curve can be placed in infinitely many one-parameter famihes; given polynomials of 
sufficiently high degree we can force any number of curves to lie in two distinct families. 
Thus it is not surprising that we run into such problems when we amalgamate. When 
we remove the repeated curves, the Pooled and Unpooled Two-Sample t-Procedures still 
give t-statistics exceeding 20 with over 200 degrees of freedom, indicating the two means 
significantly differ and supporting the claim that the repulsion decreases with increasing 
conductor 
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Figure 7: First normalized zero above the central point. 
209 rank curves from 14 rank one-parameter families, 
log(cond) S [3.26, 9.98], median = 1.35, mean = 1.36 
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Figure 8: First normalized zero above the central point. 
996 rank curves from 14 rank one-parameter families, 
log(cond) e [15.00, 16.00], median = .81, mean = .86. 
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Table 1. First normalized zero above the central point for 14 one- 
parameter families of elliptic curves of rank over Q (smaller conduc- 
tors) 



Family 


iVlcUltlll /X 


IVT^cin 1 1 
iVlcdll /i 






UlllUci 


1 


[0,1,1,1,T] 


1 78 
I .Zo 


I.:) -3 


U.ZO 


IJ.yj, V.OOJ 


7 


2 


[1,0,0,1,T] 


1 "^Q 


1 AO 


u.zy 




1 1 
1 1 


3 


[1,0,0,2,T] 


1.40 


1.41 


0.33 


[5 37 9 971 


11 


4 


[1,0,0,-1,T] 


1.50 


1.42 


0.37 


[4.70, 9.98] 


20 


5 


[1,0,0,-2,T] 


1.40 


1.48 


0.32 


[4.95, 9.85] 


11 


6 


[1,0,0,T,0] 


1.35 


1.37 


0.30 


[4.74, 9.97] 


44 


7 


[1,0,1,-2,T] 


1.25 


1.34 


0.42 


[4.04, 9.46] 


10 


8 


[1,0,2,1,T] 


1.40 


1.41 


0.33 


[5.37,9.97] 


11 


9 


[1,0,-1, 1,T] 


1.39 


1.32 


0.25 


[7.45, 9.96] 


9 


10 


[1,0,-2,1,T] 


1.34 


1.34 


0.42 


[3.26, 9.56] 


9 


11 


[1,1,-2,1,T] 


1.21 


1.19 


0.41 


[5.73, 9.92] 


6 


12 


[1,1,-3,1,T] 


1.32 


1.32 


0.32 


[5.04, 9.98] 


11 


13 


[1,-2,0,T,0] 


1.31 


1.29 


0.37 


[4.73,9.91] 


39 


14 


[-1,1,-3,1,T] 


1.45 


1.45 


0.31 


[5.76, 9.92] 


10 


All Curves 


1.35 


1.36 


0.33 


[3.26,9.98] 


209 


Distinct Curves 


1.35 


1.36 


0.33 


[3.26,9.98] 


196 



Table 2. First normalized zero above the central point for 14 one- 
parameter families of elliptic curves of rank over Q (larger conductors) 



Family 


Median 


Mean 


StDev 


log(conductor) 


Number 


1 


[0,1,1,1,T] 


0.80 


0.86 


0.23 


[15.02, 15.97] 


49 


2 


[1,0,0,1,T] 


0.91 


0.93 


0.29 


[15.00, 15.99] 


58 


3 


[1,0,0,2,T] 


0.90 


0.94 


0.30 


[15.00, 16.00] 


55 


4 


[1,0,0,-1,T] 


0.80 


0.90 


0.29 


[15.02, 16.00] 


59 


5 


[1,0,0,-2,T] 


0.75 


0.77 


0.25 


[15.04, 15.98] 


53 


6 


[1,0,0,T,0] 


0.75 


0.82 


0.27 


[15.00, 16.00] 


130 


7 


[1,0,1,-2,T] 


0.84 


0.84 


0.25 


[15.04, 15.99] 


63 


8 


[1,0,2,1,T] 


0.90 


0.94 


0.30 


[15.00, 16.00] 


55 


9 


[1,0,-1, 1,T] 


0.86 


0.89 


0.27 


[15.02, 15.98] 


57 


10 


[1,0,-2,1,T] 


0.86 


0.91 


0.30 


[15.03, 15.97] 


59 


11 


[1,1,-2,1,T] 


0.73 


0.79 


0.27 


[15.00, 16.00] 


124 


12 


[1,1,-3,1,T] 


0.98 


0.99 


0.36 


[15.01, 16.00] 


66 


13 


[1,-2,0,T,0] 


0.72 


0.76 


0.27 


[15.00, 16.00] 


120 


14 


[-1,1,-3,1,T] 


0.90 


0.91 


0.24 


[15.00, 15.99] 


48 


All Curves 


0.81 


0.86 


0.29 


[15.00,16.00] 


996 


Distinct Curves 


0.81 


0.86 


0.28 


[15.00,16.00] 


863 
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4.3.2. Rank 2 Curves. The previous results were for well-separated ranges of conductors. 
As the first set often has very small log-conductors, it is possible those values are anoma- 
lous. We therefore study two sets of curves where the log-conductors, while different, are 
close in value. The goal is to see if we can detect the effect of slight differences in the 
log-conductors on the repulsions. 

Table 3 provides the data from an analysis of 21 rank one-parameter families of elliptic 
curves over Q. The families are from IFe2l . In each family t ranges from —1000 to 1000. 
We searched for rank 2 curves with log-conductor in [15, 16]. While we study rank 2 
curves from families of rank 2 over Q in ^4.41 there the conductors are so large that we 
can only analyze a few curves in each family. In particular, there are not enough curves 
in one family with conductors approximately equal to detect how slight differences in the 
log-conductors effect the repulsions. 

Table 3 . First normalized zero above the central point for rank 2 curves 
from one-parameter families of rank over Q. The first set are curves 
with log(cond) G [15, 15.5); the second set are curves with log(coiid) £ 
[15.5, 16]. Median = pt. Mean = ji. Standard Deviation (about the Mean) 



Family 








Number 


A' 


A« 




Number 


1 


[0,1,3,1,T] 


1.59 


1.83 


0.49 


8 


1.71 


1.81 


0.40 


19 


2 


[1,0,0,1,T] 


1.84 


1.99 


0.44 


11 


1.81 


1.83 


0.43 


14 


3 


[1,0,0,2,T] 


2.05 


2.03 


0.26 


16 


2.08 


1.94 


0.48 


19 


4 


[1,0,0,-1,T] 


2.02 


1.98 


0.47 


13 


1.87 


1.94 


0.32 


10 


5 


[1,0,0,T,0] 


2.05 


2.02 


0.31 


23 


1.85 


1.99 


0.46 


23 


6 


[1,0,1,1,T] 


1.74 


1.85 


0.37 


15 


1.69 


1.77 


0.38 


23 


7 


[1,0,1,2,T] 


1.92 


1.95 


0.37 


16 


1.82 


1.81 


0.33 


14 


8 


[1,0,1,-1,T] 


1.86 


1.88 


0.34 


15 


1.79 


1.87 


0.39 


22 


9 


[1,0,1,-2,T] 


1.74 


1.74 


0.43 


14 


1.82 


1.90 


0.40 


14 


10 


[1,0,-1, 1,T] 


2.00 


2.00 


0.32 


22 


1.81 


1.94 


0.42 


18 


11 


[1,0,-2,1,T] 


1.97 


1.99 


0.39 


14 


2.17 


2.14 


0.40 


18 


12 


[1,0,-3,1,T] 


1.86 


1.88 


0.34 


15 


1.79 


1.87 


0.39 


22 


13 


[1,1,0,T,0] 


1.89 


1.88 


0.31 


20 


1.82 


1.88 


0.39 


26 


14 


[1,1,1,1,T] 


2.31 


2.21 


0.41 


16 


1.75 


1.86 


0.44 


15 


15 


[1,1,-1,1,T] 


2.02 


2.01 


0.30 


11 


1.87 


1.91 


0.32 


19 


16 


[1,1,-2,1,T] 


1.95 


1.91 


0.33 


26 


1.98 


1.97 


0.26 


18 


17 


[1,1,-3,1,T] 


1.79 


1.78 


0.25 


13 


2.00 


2.06 


0.44 


16 


18 


[1,-2,0,T,0] 


1.97 


2.05 


0.33 


24 


1.91 


1.92 


0.44 


24 


19 


[-1,1,0,1,T] 


2.11 


2.12 


0.40 


21 


1.71 


1.88 


0.43 


17 


20 


[-1,1,-2,1,T] 


1.86 


1.92 


0.28 


23 


1.95 


1.90 


0.36 


18 


21 


[-1,1,-3,1,T] 


2.07 


2.12 


0.57 


14 


1.81 


1.81 


0.41 


19 


All Curves 


1.95 


1.97 


0.37 


350 


1.85 


1.90 


0.40 


388 


Distinct Curves 


1.95 


1.97 


0.37 


335 


1.85 


1.91 


0.40 


366 



We split these rank 2 curves from the 21 one-parameter families of rank over Q into 
two sets, those curves with log-conductor in [15, 15.5) and in [15.5, 16]. We compared the 
two sets to see if we could detect the decrease in repulsion for such small changes of the 
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conductor. We have 21 families, with 350 curves in the small conductor set and 388 in the 
large conductor set. 

Remark 4.3. The families are not independent: there are 15 curves that occur twice in the 
small conductor set, and 22 in the larger. In our amalgamations of the families we consider 
both the case where we do not remove these curves, as well as the case where we do. There 
is no significant difference in the results (the only noticeable change in the table is for the 
mean for the larger conductors, which increases from 1.9034 to 1.9052 and thus is rounded 
differently). See also Remark l4!2l 

The medians and means of the small conductor set are greater than those from the large 
conductor set. For all curves the Pooled and Unpooled Two-Sample t-Procedures give t- 
statistics of 2.5 with over 600 degrees of freedom; for distinct curves the Pooled i-statistic 
is 2.16 (respectively, the Unpooled ^-statistic is 2.17) with over 600 degrees of freedom. 
As the degrees of freedom is so large, we may use the Central Limit Theorem. As there is 
about a 3% chance of observing a z-statistic of 2. 16 or greater, the results provide evidence 
against the null hypothesis (that the means are equal) at the .05 confidence level, though 
not at the .01 confidence level. 

While the data suggests the repulsion decreases with increasing conductor, it is not as 
clear as our earlier investigations (where we had ^-values greater than 10). This is, of 
course, due to the closeness of the two ranges of conductors. We apply non-parametric 
tests to further support our claim that the repulsion decreases with increasing conductors. 
For each family in Table 3, write a plus sign if the small conductor set has a greater mean 
and a minus sign if not. There are four minus signs and seventeen plus signs. The null 
hypothesis is that each mean is equally likely to be larger. Under the null hypothesis, the 
number of minus signs is a random variable from a binomial distribution with = 21 and 
6 = ^. The probability of observing four or fewer minus signs is about 3.6%, supporting 
the claim of decreasing repulsion with increasing conductor For the medians there are 
seven minus signs out of twenty-one; the probability of seven or fewer minus signs is 
about 9.4%. Every time the smaller conductor set had the lesser mean, it also had the 
lesser median; the mean and median tests are not independent. 
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4.4. One-Parameter Families of Rank 2 Over Q. 

4.4.1. Family = — T^x + T^. We study the first normalized zero above the central 
point for 69 rank 2 elliptic curves from the one-parameter family = x^ — T^x + 
of rank 2 over Q. There are 35 curves with log(cond) e [7.8, 16.1] in Figure 9 and 34 
with log(cond) G [16.2, 23.3] in Figure 10. Unlike the previous examples where we chose 
many curves of the same rank from different families, here we have just one family. As the 
conductors grow rapidly, we have far fewer data points, and the range of the log-conductors 
is much greater However, even for such a small sample, the repulsion decreases with in- 
creasing conductors, and the shape begins to approach the conjectured distribution. The 
Pooled and Unpooled Two-Sample t-Procedures give t-statistics over 5 with over 60 de- 
grees of freedom, and we may use the Central Limit Theorem. As the probability of a 
z-value of 5 or more is less than lO^'* percent, the data does not support the null hypoth- 
esis (i.e., the data supports our conjecture that the repulsion decreases as the conductors 
increase). 
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Figure 9: First normalized zero above the central point 
from rank 2 curves in the family y^ = x^ — T'^x + T^. 
35 curves, log(cond) G [7.8, 16.1], median = 1.85, mean = 1.92, 
standard deviation about the mean = .41 
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Figure 10: First normalized zero above the central point 
from rank 2 curves in the family y"^ ^ x^ — T^x + T^. 
34 curves, log(coiid) e [16.2, 23.3], median = 1.37, mean = 1.47, 
standard deviation about the mean = .34 
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4.4.2. Rank 2 Curves. We consider 21 one-parameter families of rank 2 over Q, and in- 
vestigate curves of rank 2 in these families. The families are from IFe2l . We again amal- 
gamated the different families, and summarize the results in Table 4. 

The difference between these experiments and those of ^4. 3. 21 is that, while both deal 
with one-parameter families over Q, here we study curves of rank 2 from families of rank 
2 over Q; earlier we studied curves of rank 2 from families of rank over Q. If the Density 
Conjecture (with orthogonal groups) is correct for the entire one-parameter family, in the 
limit 0% of the curves in a family of rank r have rank r + 2 or greater Thus our previous 
investigations of curves of rank 2 in a family of rank over Q were a study of a measure 
zero subset. Unlike curves of rank 2 in families of rank 2 over Q, we have no theoretical 
evidence supporting a proposed random matrix model for curves of rank 2 in families of 
rank 0. We compare the results from rank 2 curves in rank 2 families over Q to the rank 2 
curves from rank families over Q in ^.51 

Table 4. First normalized zero above the central point for 21 one- 
parameter families of rank 2 over Q with log(cond) G [15, 16] and 
t G [0, 120]. The median of the first normalized zero of the 64 curves is 
1.64. 



Family 


Mean 


Standard Deviation 


log(conductor) 


Number 


1 


[l,T,0,-3-2T,l] 


1.91 


0.25 


[15.74,16.00] 


2 


2 


[1,T,-19,-T-1,0] 


1.57 


0.36 


[15.17,15.63] 


4 


3 


[1,T,2,-T-1,0] 


1.29 




[15.47, 15.47] 


1 


4 


[1,T,-16,-T-1,0] 


1.75 


0.19 


[15.07,15.86] 


4 


5 


[1,T,13,-T-1,0] 


1.53 


0.25 


[15.08,15.91] 


3 


6 


[1,T,-14,-T-1,0] 


1.69 


0.32 


[15.06,15.22] 


3 


7 


[1,T,10,-T-1,0] 


1.62 


0.28 


[15.70,15.89] 


3 


8 


[0,T,11,-T-1,0] 


1.98 




[15.87,15.87] 


1 


9 


[1,T,-11,-T-1,0] 










10 


[0,T,7,-T-1,0] 


1.54 


0.17 


[15.08,15.90] 


7 


11 


[l,T,-8,-T-l,0] 


1.58 


0.18 


[15.23,25.95] 


6 


12 


[1,T,19,-T-1,0] 










13 


[0,T,3,-T-1,0] 


1.96 


0.25 


[15.23, 15.66] 


3 


14 


[0,T,19,-T-1,0] 










15 


[1,T,17,-T-1,0] 


1.64 


0.23 


[15.09, 15.98] 


4 


16 


[0,T,9,-T-1,0] 


1.59 


0.29 


[15.01, 15.85] 


5 


17 


[0,T,1,-T-1,0] 


1.51 




[15.99, 15.99] 


1 


18 


[l,T,-7,-T-l,0] 


1.45 


0.23 


[15.14, 15.43] 


4 


19 


[1,T,8,-T-1,0] 


1.53 


0.24 


[15.02, 15.89] 


10 


20 


[l,T,-2,-T-l,0] 


1.60 




[15.98, 15.98] 


1 


21 


[0,T,13,-T-1,0] 


1.67 


0.01 


[15.01, 15.92] 


2 


All Curves 


1.61 


0.25 


[15.01, 16.00] 


64 



Remark 4.4. There are 23 rank 4 curves in the 21 one-parameter families of rank 2 over 
Q with log-conductors in [15, 16] and t G [0, 120]. For the first normalized zero above the 
central point, the median is 3.03, the mean is 3.05, and the standard deviation about the 
mean is 0.30. 
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4.5. Comparison Between One-Parameter Families of Different Rank. In Table 5 we 
investigate how the first normalized zero above the central point of rank 2 curves depends 
on how the curves are obtained. The first family is rank 2 curves from the 21 one -parameter 
famiUes of rank over Q from Table 3, while the second is rank 2 curves from the 21 
one-parameter families of rank 2 over Q from Table 4; in both sets the log-conductors 
are in [15, 16]. A t-Test on the two means gives a t-statistic of 6.60, indicating the two 
means differ. Thus the mean of the first normalized zero above the central point of rank 
2 curves in a one-parameter family over Q (for conductors in this range) depends on how 
we choose the curves. For the range of conductors studied, rank 2 curves from rank 
one-parameter families over Q do not behave the same as rank 2 curves from rank 2 one- 
parameter families over Q. 

Table 5 . First normalized zero above the central point. The first family 
is the 701 rank 2 curves from the 21 one-parameter families of rank 
over Q from Table 3 with log(cond) G [15, 16]; the second family is the 
64 rank 2 curves from the 21 one-parameter families of rank 2 over Q 
from Table 4 with log(cond) g [15, 16]. 



Family 


Median 


Mean 


Std. Dev. 


Number 


Rank 2 Curves, Families Rank over Q 
Rank 2 Curves, Families Rank 2 over Q 


1.926 
1.642 


1.936 
1.610 


0.388 
0.247 


701 
64 
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4.6. Spacings between normalized zeros. For finite conductors, even when there are no 
zeros at the central point, the first normaUzed zero above the central point is repelled from 
the predicted N oo scaling limits. The repulsion increases with the number of zeros 
at the central point and decreases with increasing conductor For an elliptic curve E, let 
zi, 22, Z3, . . . denote the imaginary parts of the normalized zeros above the central point. 
We investigate whether or not Zj+i — Zj depends on the repulsion from the central point. 
We consider the following two sets of curves in Table 6: 

• the 863 distinct rank curves with log(cond) G [15, 16] from the 14 one-parameter 
families of rank over Q from Table 2; 

• the 701 distinct rank 2 curves with log(cond) G [15, 16] from the 21 one-parameter 
families of rank over Q from Table 3. 

In Table 6 we calculate the median and mean for Z2 — zi, — z^ and Z3 — zi. The last 
statistic involves the sum of differences between two adjacent normalized zeros, and allows 
the possibility of some effects being averaged out. While the normalized zeros are repelled 
from the central point (and by different amounts for the two sets), the differences between 
the normalized zeros are statistically independent of this repulsion. We performed a i-Test 
on the means in the three cases. For each case the f-statistic was less than 2, strongly 
supporting the null hypothesis that the differences are independent of the repulsion. 

Table 6. Spacings between normalized zeros. All curves have 
log(cond) G [15, 16], and Zj is the imaginary part of the j* normal- 
ized zero above the central point. The 863 rank curves are from the 
14 one-parameter families of rank over Q from Table 2; the 701 rank 
2 curves are from the 21 one-parameter families of rank over Q from 
Table 3. 





863 Rank Curves 


701 Rank 2 Curves 


t-Statistic 


Median zi 


- 21 


1.28 


1.30 




Mean zi 


- 21 


1.30 


1.34 


-1.60 


StDev Z2 


- 


0.49 


0.51 




Median Z3 


~ Z2 


1.22 


1.19 




Mean Z3 


- Z2 


1.24 


1.22 


0.80 


StDev Z3 


- Z2 


0.52 


0.47 




Median Z3 


- Zl 


2.54 


2.56 




Mean Z3 


- Zl 


2.55 


2.56 


-0.38 


StDev Z3 


- Zl 


0.52 


0.52 





We have consistently observed that the more zeros at the central point, the greater the 
repulsion. One possible explanation is as follows: for rank 2 curves in a rank one- 
parameter family over Q, the first zero above the central point collapses down to the central 
point, and the other zeros are left alone. As the zeros are symmetric about the central point, 
the effect of one zero above the central point collapsing is to increase the number of zeros 
at the central point by 2. 

For our 14 one-parameter families of elliptic curves of rank over Q and log-conductors 
in [15, 16], we studied the second and third normalized zero above the central point. The 
mean of the second normalized zero is 2.16 with a standard deviation of .39, while the 
third normalized zero has a mean of 3.41 and a standard deviation of .41. These numbers 
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Statistically differ from the first and second normalized zeros of the rank 2 curves from 
our 21 one-parameter families of rank over Q with log-conductor in [15, 16], where the 
means were respectively 1.93 (with a standard deviation of .39) and 3.27 (with a standard 
deviation of .39). Thus while for a given range of log-conductors the average second 
normalized zero of a rank curve is close to the average first normalized zero of a rank 2 
curve, they are not equal and the additional repulsion from extra zeros at the central point 
cannot be entirely explained by only collapsing the first zero to the central point while 
leaving the other zeros alone. 

Remark 4.5. As the second (resp., third) normalized zero for rank curves in rank 
families over Q is 2. 16 (resp., 3.41) while the first (resp., second) normalized zero for rank 
2 curves in rank families over Q is 1.93 (resp., 3.27), one can interpret the effect of the 
additional zeros at the central point as an attraction. Specifically, for curves of rank 2 
in a rank family over Q, by symmetry two zeros collapse to the central point, and the 
remaining zeros are then attracted to the central point, being closer than the corresponding 
zeros from rank curves. As remarked in §3.5 of tFari), the term "lowest zero" is not well 
defined when there are multiple zeros at the central point. We can either mean the first zero 
above the central point, or one of the many zeros at the central point. In all cases, for finite 
conductors there is repulsion from the N ^ oo scaling limits of random matrix theory; 
however, "attraction" might be a better term for the effect of additional zeros at the central 
point, though the current terminology is to talk about repulsion of zeros at the central point. 

We now study the differences between normalized zeros coming from one-parameter 
families of rank 2 over Q. Table 7 shows that while the normalized zeros are repelled from 
the central point, the differences between the normalized zeros are statistically independent 
of the repulsion. We performed a i-Test for the means in the three cases studied. For two of 
the three cases the i-statistic was less than 2 (and in the third it was only 2.05), supporting 
the null hypothesis that the differences are independent of the repulsion. 

Table 7. Spacings between normalized zeros. All curves have 
log(cond) G [15, 16], and Zj is the imaginary part of the normal- 
ized zero above the central point. The 64 rank 2 curves are the 21 one- 
parameter families of rank 2 over Q from Table 4; the 23 rank 4 curves 
are the 21 one-parameter families of rank 2 over Q from Table 4. 





64 Rank 2 Curves 


23 Rank 4 Curves 


t-Statistic 


Median Z2 


- Zl 


1.26 


1.27 




Mean Z2 


- 21 


1.36 


1.29 


0.59 


StDev Z2 


- Zl 


0.50 


0.42 




Median zs 


- Z2 


1.22 


1.08 




Mean z^, 


- Z2 


1.29 


1.14 


1.35 


StDev 23 


- Z2 


0.49 


0.35 




Median 23 


- Zl 


2.66 


2.46 




Mean 23 


- Zl 


2.65 


2.43 


2.05 


StDev Z3 


- Zl 


0.44 


0.42 





We performed one last experiment on the differences between normalized zeros. In 
Table 8 we compare two sets of rank 2 curves: the first are the 21 one-parameter families 

^^The Pooled and Unpooled t-statistics in both experiments are greater than 6, providing evidence against the 
null hypothesis that the two means are equal. 
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Table 8. Spacings between normalized zeros. All curves have 
log(cond) G [15, 16], and Zj is the imaginary part of the j* normal- 
ized zero above the central point. The 701 rank 2 curves are the 21 
one-parameter families of rank over Q from Table 3, and the 64 rank 2 
curves are the 21 one-parameter families of rank 2 over (Q) from Table 4. 





701 Rank 2 Curves 


64 Rank 2 Curves 


t-Statistic 


Median Z2 


- Zl 


1.30 


1.26 




Mean Z2 


- Zl 


1.34 


1.36 


0.69 


StDev 22 


- Zl 


0.51 


0.50 




Median z^ 


- Z2 


1.19 


1.22 




Mean 23 


- Z2 


1.22 


1.29 


1.39 


StDev 23 


- Z2 


0.47 


0.49 




Median 23 


- Zl 


2.56 


2.66 




Mean 23 


- Zl 


2.56 


2.65 


1.93 


StDev 23 


- Zl 


0.52 


0.44 





of rank over Q from Table 3, while the second are the 21 one-parameter famiUes of rank 2 
over Q from Table 4. While the first normalized zero is repelled differently in the two cases, 
the differences are statistically independent from the nature of the zeros at the central point, 
as indicated by all i-statistics being less than 2. This suggests that the spacings between 
adjacent normalized zeros above the central point is independent of the repulsion at the 
central point; in particular, this quantity does not depend on how we construct our family 
of rank 2 curves. 

5. Summary and future work 

As the conductors tend to infinity, theoretical results support the validity of the N ^ 00 
scaling limit of the Independent Model for all curves in one-parameter families of ellip- 
tic curves of rank r over Q; however, it is unknown what the correct model is for the 
sub-family of curves of rank r + 2. The experimental evidence suggests that the first 
normalized zero, for small and finite conductors, is repelled by zeros at the central point. 
Further, the more zeros at the central point, the greater the repulsion; however, the repul- 
sion decreases as the conductors increase, and the difference between adjacent normalized 
zeros is statistically independent of the repulsion and the rank of the curves. 

At present we can calculate the first normalized zero for log-conductors about 25. While 
we can use more powerful computers to study larger conductors, it is unlikely these con- 
ductors will be large enough to see the predicted limiting behavior It is interesting that, 
unlike the excess rank investigations, we see noticeable convergence to the limiting theo- 
retical results as we increase the conductors. 

An interesting project is to determine a theoretical model to explain the behavior for 
finite conductors. In the large-conductor limit, analogies with the function field and cal- 
culations with the explicit formula lead us to the Independent Model for curves of rank r 
from families of rank r over Q, and theoretical results in the number field case support this. 
It is not unreasonable to posit that in the finite-conductor analogues the size of the matrices 
should be a function of the log-conductors. Unfortunately the statistics for the finite N xN 
random matrix ensembles are expressed in terms of eigenvalues of integral equations, and 
are usually only plotted in the N 00 scaling limit. This makes comparison with the 
experimental data difficult, and a future project is to analyze the finite N cases by using 
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the finite N kernels. Such an analysis will facilitate comparing the finite N limits of the 
Independent and Interaction Models for curves of rank r + 2 from families of rank r over 

Q. 
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Appendix A. "Harder" Ensembles of Orthogonal Matrices 
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In this appendix we derive the conditional (interaction) eigenvalue probability mea- 
sure (I2.13> and illustrate how it affects eigenvalue statistics near the central point 1, in 
particular through repulsion (observed via the 1-level density). We also explain the re- 
lation to the classical Bessel kernels of random matrix theory, and to other central-point 
statistics. 

A. 1 . Full Orthogonal Ensembles. In view of our intended application we will be con- 
cerned exclusively with random matrix ensembles of orthogonal matrices in what follows. 
If we write the eigenvalues (in no particular order) of a special'^ orthogonal matrix of 
size 2N (resp., 2N + 1) as {ie'^^jf (resp., {+1} U {±e*''^ }f ) with Q < Oj < n then 
the iV-tuple Q = . . . , 0^) parametrizes the eigenvalues. In terms of the angles 9j, the 
probability measure of the eigenvalues induced from normalized Haar measure on SO(2 A^) 
(resp., on SO(2A^ + 1) upon discarding one forced eigenvalue of +1) can be identified with 
a measure on [0, tt]^, 

deo(e) = C'^^ \Y (cos ^fc - cos 6*^)2 J| dO^ (A.22) 

i<i<fc<JV i<i<JV 

dei{<d) = C^^^ {cosOk-cosOjf J| sm^{%)dej (A.23) 

i<i<fc<JV i<i<JV 

in the 2A^ and 2A^ + 1 cases, respectively, as shown in IConllKaSal l: the normalization 
constants C'^'' ensure that the measures on the right-hand sides are probability measures. 
Note that formulas ( IA.22> and (IA.23> are symmetric upon permuting the 9j\, so issues 
related to a choice of a particular ordering of the eigenvalues of the matrix are irrelevant. 
More importantly, observe the quadratic exponent of the differences of the cosines. 

The statistical behavior of the eigenvalues near +1 is closely related to the order of 
vanishing of the measures above at = 0. We change variables and replace the eigenvalues 
g±(6i j-jjg levels 

X = cos 6* (A.24) 
so the measures above become measures on [—1, +1]^: 

N 

deoiX) = J] {xk - x,)^ x,)-i{l + x,)-idx, (A.25) 

i<j<k<N j=l 
N 

dei{X) = Cl^'^l[ixk-x,fl[{l-x^)i(l + x,)-idx,, (A.26) 

3<k 3 = 1 



'That is, of determinant one. 
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where X = {xi , . . . , xn) and C^'"* are suitable normalization constants. Here we observe 
the appearance of the weight functions on [—1, 1] 

.(.) = (l-.r(l + .r^ a=h/' (A.27) 
^ ^ ^ ^ ' I +1/2 for SO(2A^ + 1). 



By the Gaudin-Mehta theory (see for example IMehl '). and in view of the quadratic expo- 
nent of the differences of the "levels" Xj, the study of eigenvalue statistics using classical 
methods is intimately related to the sequence of orthogonal polynomials with respect to 
the weight w{x)}^ 

In classical random matrix theory terminology (especially in the context of the Laguerre 
and Jacobi ensembles) the endpoints —1, +1 are called the "hard edges" of the spectrum 
because the probability measure, considered on R^, vanishes outside [—1, +1]^. We will 
keep calling = 0, tt the "central points" (endpoints of the diameter with respect to which 
the spectrum is symmetric). Phenomena about central points and hard edges are equivalent 
in view of the change of variables ( IA.24> . Perhaps not surprisingly, the parameter a, which 
dictates the order of vanishing of the weight function 'w{x) at the hard edge +1, suffices 
to characterize the mutually different statistics near the central point in each of SO(even) 
and SO(odd). However, the importance of this parameter is best understood in the context 
of certain sub-ensembles of SO as described below. 

A.2. Conditional ("Harder") Ortliogonal Ensembles. The conditional eigenvalue mea- 
sure for the sub-ensemble S'O'^^'') (2^) of SO{2N) consisting of matrices for which some 
2r of the 2N eigenvalues are equal to +1 can easily be obtained from (IA.25t . Let 

n i^k~x,f n w{x,) (A.28) 

l<j<k<N i<j<N 

be the normalized probability density function of the levels for S0{2N), where w{x) is as 
in ( IA.27> with a = —1/2 and m ~ 0. Now let ti, ... ,tr be chosen so < < 1, let 
K = Yiji^ ^ 1] ^ii'l I = J X K for some box J C [—1, l]^"*". This means we are 
constraining r pairs of levels to lie in a neighborhood of a:; = 1 (or equivalently that we 
are construing r pairs of eigenvalues to lie in circular sectors about the point 1 on the unit 
circle). Thus, the conditional probability that the remaining N — r pairs of eigenvalues lie 
in J is given by 

f{x)dx 

f{x)dx 
'[-i,i]"-''xa: ' 

where T ~ (ti, . . . , tr). The conditional probabiUty measure of the eigenvalues for the 
sub-ensemble S'0(20(2iV) is the limit as all tk 0+ of F{t; J), as a function of the box 
J, call it G{J). Applying L'Hopital's rule r times to the quotient (IA.29> (once on each 
variable tk) and using the fundamental theorem of calculus we get 

, . SA^MX))\M{X,T)fw{X)dX ■ (Van(r))^^(T) 

T™+/j^^^j„_(Van(X))2(Af(X,r))2u;(X)dX.(Van(T))2u;(r) ^ ' ' 



F{T;J) = ''■^''^ , (A.29) 



'■'The inner product being (/, 3) = f(x)g(x)'w(x)dx. 
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where X = {xi, . . . , XN~r), and 

Van(X)= Yl (xk-x,) Van(T)= [] {tk - tj) 

l<j<k<N~r i<j<k<r 

w{x)= n ^(^^o n ^(i-^fc) 

l<j<N-r l<k<r 



\<j<N -r 
l<k<T 



Naturally, the factors of Van(r) and cancel in equation ( IA.30> . Since M{X, T) 

is bounded, the integrands in equation ( lA.BOt are uniformly dominated by an integrable 



function, ensuring that we can let all tj — > in the integrands of ( IA.30> to obtain 

_ //Van(X))^(M(X,0))^^(X)rfX 
^ ' /[o_i]„_(Van(X))2(M(X,0))2ti;(X)dX- ^ ' ' 

Now observe that {M{X,{))Yw{X) = ni<j<r^(^i) where 'w{x) is given by equa- 
tion JA.27> with a replaced by a = a + 2, so the probability measure of the eigenvalues 
for S'O'-^'') (2iV) is obtained from that of SO{2{N — ra)) simply by changing the weight 

function w ^ w. Explicitly, the probability measure of the eigenvalues of the ensemble 

S'Cl(2r)(2A^) is given 

ds^iX) = d;pl^\[{xk - x,f\{{l - z,)—^ (A.32) 

j<k j j 

where m = 2r, X = {xi , . . . , XN-r), the indices j, k range from 1 to — r, and C^}^ 
are suitable normalization constants (equal to the reciprocal of the denominator of the 



right-hand side of equation ( IA.31> .) 

The same argument shows that the sub-ensemble SO'^^''+^^2N + 1) of S0{2N + 1) 
consisting of matrices for which 2r + 1 eigenvalues are equal to + 1 has the same eigenvalue 
measure with m = 2r + 1. Because the density of the measure vanishes to a higher 

order near the edge +1 the larger m is, we will say that the edge becomes harder when m 
is larger (whence the title of this section), and call m its hardness. 



A. 3. Independent Model. It is important to observe that the presence of the m-multiple 
eigenvalues at the central point in these harder sub-ensembles of orthogonal matrices has 
a strong repelling effect due to the extra factor (1 — a;)™ multiplied by the weight (1 — 
2 (1 + x)^3 of SO(even). For comparison purposes consider the following situation, 
first in the SO(even) case. The number of eigenvalues equal to +1 of any S0{2N) matrix 
is always an even number 2r, and one may consider the sub-ensemble A2N,2r of SO{2N), 

A2N,2r = I (^^'''""' J : 3 € SO{2N - 2r)| , (A.33) 

which is just S0{2N — 2r) in disguise. This is certainly a sub-ensemble of SO{2N) 
consisting of matrices with at least 2r eigenvalues equal to +1, albeit quite a different one 
from the 2r-hard sub-ensemble of SO{2N) described before. For example, the eigenvalue 
measure (apart from the point masses at the last 2r eigenvalues) for A2N,2r is 



deo{xi, . . . ,XN-r), 



(A.34) 
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and not de2r{xi, . . . , XN^r)- The same observation applies in the S0{2N +1) case: with 
the obvious notation, the eigenvalue measure for A2N+i,2r is'"* 

(iei(0i,...,0w_,), (A.35) 

and not de2r+i(^i, ■ ■ ■ , On-t)- 



A.4. 1 -Level Density: Full Orthogonal. Before dealing with the harder sub-ensembles, 
we make some comments about the hard edges of the full SO(even) and SO(odd). The 
local statistics near the point +1 are dictated by the even '+' (resp., odd ' — ') Sine Kernels 



^±(C,^/)-^(C,'/)±5(e,-r7) (A.36) 

in the case of SO(even) (resp., SO(odd)); see IIKaSall fKaSa2l . Here ^,?7 are rescaled 
variables centered about the value 0, namely related to the original variables by'^ 

X = cos (^^) (A.37) 

and S{x, y) = sin(7ra;) / (ttx) is the Sine Kernel, which has the universal property of de- 
scribing the local statistics at any bulk point of any ensemble with local quadratic local 
level repulsion IIDKMcVZI . However, it is not the Sine Kernel but its even (resp., odd) 
counterparts that dictate the local statistics near the central point. For example, the central 
one-level density is given by the diagonal values at a; = j/ of the respective kernel: 

sm Sttx 

p+{x) = 1 + — , for SO(even), (A.38) 

siu Sttx 

p-{x) ^ 1 + — + S(x), for SO(odd). (A.39) 

znx 

(In the SO(odd) case the Dirac delta reflects the fact that any such matrix has an eigenvalue 
at the central point.) Observe that /?_ vanishes to second order, whereas does not vanish 
at the central point x = 0. 



A. 5. 1-Level Density: Harder Orthogonal. We return to the more general case of m- 
hard ensembles of orthogonal matrices. Because the classical Jacobi polynomials {pjf''^ }q 
are orthogonal with respect to the weight w{x) = (1 — a;)''(l + x)^ on [—1, 1], the local 
statistics near the central point x ^ +1 are derived from the asymptotic behavior of these 
polynomials at the right edge of the interval [— 1 , + 1] . More specifically, the relevant ker- 
nel which takes the place of the (even or odd) Sine Kernel is the "edge limit" as ^ oo of 
the Christoffel-Darboux/Szego projection kernel K^^'^'^ {x, y) onto polynomials of degree 
less than N in 1, 1], (1 — xY{l + x)^dx) (via the change of variables ( IA.37» . For 



Observe that a matrix in A-2N+i,2r has 2r + 1 eigenvalues equal to +1 and N — r other pairs eigenvalues. 

'^This is justified by the fact N/ir is the average (angular) asymptotic density of the eigen-angles Bj of a 
random orthogonal matrix, hence asymptotic equidistribution holds — away from the central points ! 

""If the central point were not atypical, then the local density would be dictated by the diagonal values 
S{x,x) = 1 of the Sine Kernel. 

'^This "edge limit" and the ensuing Bessel Kernels are also observed in the somewhat simpler context of the 
so-called (unitary) Laguen'e ensemble. 
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Figure 1. The central 1-level densities (dotted) and p_ (dash- 
dotted) versus the "bulk" 1-level density p = 1 observed away from 
the central points. 

the edge +1, the limit depends only on the parameter a and is equal to the Bessel kernel"^ 

^'"^(^'^) = ^F^[<'^a+lK)Ja(7rr,) - JaiTTOTTVJa+liTTv)], (A.40) 
B^^H^^O = f K)[^a«) - ^a-l«)Ja+lK)], (A.41) 

where J^, stands for the Bessel function of the first kind INWiO . 

It is a little more natural for our purposes to use the hardness to, rather than a = m — ^, 
as the parameter, so we define 

(A.42) 

k^'^Hx^y) = -K'-"'Hx/T^,y/T^)- (A-43) 

TT 

Using the recursion formula for Bessel functions we obtain an alternate formula to ( IA.41> 
for the diagonal values of the kernel: 

/c('")(x,x) ^ |[J,„+i + J„,_i{x)'] ~ (to- i)J,„+i (x)J,„_|(x). (A.44) 

Except for to times a point mass at x = 0, the central m-hard 1-level density is given by 

p„,{x)=K^"'Hx,x). (A.45) 



A. 6. Spacing Measures. In this section we state some well-known formulas giving the 
spacing measures or "gap probabilities" at the central point. Their derivation is stan- 
dard and depends only on knowledge of the edge limiting kernels A'*^'"' (see for instance 
IKaSalllMehllTWl ). Let s) be the Umit, as iV cx), of the probabihty t hat ex - 

actly k of the ^j's lie on the interval (0, s), where £_j is related to xj via equation ( IA.37t . 
Also let p^™\k; s)ds be the conditional probability that the (fc + l)-st of the ^^'s, to the 
right of ^ = 0, lies in the interval [s, s + ds), in the limit N ^ oo. 

"^In fact, the even Sine Kernel 5"+ = _b'~ 2' whereas the odd Sine Kernel S— = B^^\ 
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Figure 2. The m-hard 1-level edge densities for m = 0, 1, . . . , 5. 



Abusing notation, let K^"^^\s denote the integral operator on i^([0, s], dx) with kernel 



If / denotes the identity operator, then the following formulas hold: 



T=-l 



(A.46) 

(A.47) 
(A.48) 



On the right-hand side of (IA.47t . 'det' is the Fredhobn determinant: for an operator with 
kernel /C, 

dct(/ + /C) = 1 + y ^ / / dct(IC(x,,Xk))dx,,...dxi. (A.49) 

^ — ' n' / nxn 

n=l 

Identical formulas hold even for finite N provided that the limiting kernel i^t™) is re- 
placed by the Christoffel-Darboux/Szego projection kernel K^_^' ^ associated to the 
weight 'w{x) of (IA.27> with a = m — ^, acting on 1, l\,w{x)dx). In this case, the 

corresponding operator is of finite rank, the Fredholm determinant agrees with the usual 
determinant, and the series \AA9l is finite. 

A. 7. Explicit Kernels. In view of the relation between Bessel Functions of the first kind 
of half-integral parameter and trigonometric functions, it is possible to write the kernels 
j<[{m) jjj tgj-jjjs of elementary functions. 



A.7.1. m ~ 0." The Even Sine Kernel. 



Ko{x,y) = S+{x,y) 



simr{x — y) sin7r(a; + y) 



7t{x - y) 



7r(a; + y) 



(A.50) 
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The one-level density is 



p+(x) = 5+(a:,a:) = l + ^|^. (A.51) 



The Fourier transform of the one-level density is 

p+{u)^S{u) + ^I{u), (A.52) 
where I{u) is the characteristic function of the interval [—1,1]. 

A. 7. 2. m = I: The Odd Sine Kernel. 

I \ Q t \ siii7r(a;-?/) sin7r(x + y) 

Ki{x,y) = S_{x,y)^ -, ^ -, — \ — ^ ■ (A.53) 

7r(x — y) tt[x + y) 



The one-level density is 



p^{x) = S_ix,x) =d{x) + l-^^^^^^. (A.54) 

ZTTX 



The Fourier transform of the one-level density is 

p^{u)^Siu) + l-^Iiu). (A.55) 

A. 7. 3. m = 2: The "Doubly Hard" Kernel. 

sm.-K{x-y) sm.Tr{x + y) sin ttx sin Try c<^^ 

K2{x,y)^ — — + — - — ■ — r 2 . (A.56) 

7r(x- — y) 7r(x + y) irx iry 



sin27ra; „ / sinTrx^ ^ 



(A.59) 



The one-level density is 

p^{x) = 2S{x) + 1 + --^^^-^ - 2 { ^^^^ ] . (A.57) 

ZTTX \ TTX / 

The Fourier transform of the one-level density is 

P2iu) = 5{u) + 2 + (2|u| - |)/(u). (A.58) 
A. 7.4. m = 3." The "Triply Hard" Kernel. 

K:i{x,y) = Ki{x,y) + [l + ^— ] Ko{x,y) 

cos TTCC cos Try sin ttx sin Try 
■nx ny (nx)'^ {t^v)'^ 
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